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Abstract 

An explicit evaluation of the two pion transverse momentum dependent parton distributions at 
leading twist is presented, in the framework of the Nambu-Jona Lasinio model with Pauli-Villars 
regularization. The transverse momentum dependence of the obtained distributions is generated 
solely by the dynamics of the model. Using these results, the so called generalized Boer-Mulders 
shift is studied and compared with recent lattice data. The obtained agreement is very encouraging, 
in particular because no additional parameter has been introduced. A more conclusive comparison 
would require a precise knowledge of the QCD evolution of the transverse momentum dependent 
parton distributions under scrutiny. 
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1 Introduction 

The three-dimensional (3D) hadronic structure in momentum space can be accessed through 
the transverse momentum dependent parton distributions (TMDs) [1], measured mainly in semi- 
inclusive deep inelastic scattering (SIDIS) or in Drell-Yan (DY) processes. For the nucleon target 
a large amount of theoretical work is being done, driven by recent and forthcoming impressive 
experimental efforts (see, e.g, laiiii] and references therein). In this paper we discuss pion TMDs, 
which are experimentally probed through the DY process (see, e.g., [5] for a recent report). 

At leading-twist, the pion structure is described in terms of two TMDs, the unpolarized one, 
k'^), describing the number density of partons with longitudinal momentum fraction x 
and transverse momentum kx, and the Boer-Mulders TMD, also called Boer-Mulders function, 
k'^ ) mni- The latter is not a density, being generated by spin-orbit correlations of trans¬ 
versely polarized partons; it is chiral-odd and therefore not accessible in DIS, and it is “naively” 
time-reversal odd, i.e., under time reversal the correlation flips the sign. 

TMDs are non perturbative quantities and they have not been calculated from first principles, 
although recently lattice data have been produced for the pion. In particular, the lattice calculation 
in Ref. [8], performed at the pion mass ttItt = 518 MeV, is the update of preliminary results reported 
in laiiQ]. Pion TMDs have been estimated also in models of the pion structure, such as spectator 
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models [niiiaiis], bag models [14], covariant model of the pion with Pauli-Villars regulators, in 
the unpolarized case |15j . and in a light-front constituent quark model m- 

In this paper we present the calculation of and in the model of Nambu 

and Jona-Lasinio (NJL) [17j . 

The NJL model is the most realistic model for the pseudoscalar mesons based on a local quantum 
field theory built with quarks. It respects the realization of chiral symmetry and gives a good 
description of low energy properties. Mesons are described as bound states, in a fully covariant 
way, using the Bethe-Salpeter amplitude, in a field theoretical framework. In this way, the Lorentz 
covariance of the problem is preserved. The NJL model is a non-renormalizable field theory and 
therefore a cut-off procedure has to be implemented. Here, the Pauli-Villars regularization scheme 
has been chosen, because it respects the gauge symmetry of the problem. The NJL model, together 
with its regularization procedure, can be regarded as an effective theory of QCD. 

The NJL model has a long tradition of successful predictions of different observables related to 
the parton structure of pseudoscalar mesons, such as the parton distribution [T8l[T9], generalized 
parton distributions [20] , distribution amplitudes [21] , transition distribution amplitudes [22l [23] , 
transition form factors [Ml 113 US]- Here, for the first time, we apply the same scheme to the 
calculation of the pion TMDs. This will permit to obtain a dynamical dependence, at variance 
with various other model analyses where its analytical trend was assumed, and to compare it with 
very recent lattice data [5] . 

The paper is structured as follows. In Section 2 we describe our approach obtaining the formal 
results. In the third Section we discuss the numerical results and, at the end, we perform the 
comparison with lattice data. Conclusions are eventually presented in the last section. 

2 TMDs in the NJL model 

For a spinless particle, such as the pion, only two leading twist TMDs arise, in contrast to the eight 
found for spin-i particles [3]. The TMD /i ^ is simply the unpolarized quark distribution, whereas 
the Boer-Mulders (BM) function [5], describes the distribution of transversely polarized quarks 
in the pion. The BM function is odd under time reversal (T-odd). A non-zero value for this function 
is originated by the final and initial state interactions, in the SIDIS and DY processes, respectively, 
which break the symmetry of the events under time reversal. 

The calculation of /i,,r and in the NJL model will be described in the following two sub¬ 
sections, respectively. 

2.1 Unpolarized TMD 

The unpolarized quark TMD in the pion is defined as follows 


A., - 2 7 (2,)3 « 

X (oo,u) CT3) £0(00,0) V'(o)b) , (I) 

where c = I(—I) stands for the u{d) case, = xP^, U = 0, and the gauge link is given by 


£ 




( 00 , C ) =£> exp (^-igs J A+ (j] ,^t) drj ^ 


( 2 ) 


with gs the strong SU (3)^ coupling constant. 

To fix the ideas, we consider the u quark TMD in a tt’*'. At zero order in gs, one gets 


/i.TT (a;, fcr) = ^ J ^ (pI V'(C ,U) 7"^ ^ (1 +'^3) V'(0) b) ■ (3) 
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Figure 1: The diagrams to be considered in the evaluation of the u TMD in a tt"*", 


Eq. ([3]). 


The two diagrams contributing to this quantity are shown in Fig. 1. The contribution of the 
diagram in the left panel is 


fl,Tr {x, krp^ 






X Tr (^q - ^ 

X -i+iSp (^q+ 


'^9'^qq r— iSp 


igTrqq *75 *"+ 



(4) 


where Tr implies traces in color, flavor and Dirac matrices, Sp{p) is the Feynman propagator and 
T± = ± iT2). The other diagram of Fig. 1, corresponding to the propagation of a cr particle, 

which gives sometimes important contributions (see for instance the calculation of pion GPDs in 
EQ]), vanishes in this case, where a diagonal matrix element of a bi-local current is involved. 

After Pauli Villars regularization (see the Appendix for details), the final result is 


/i.TT {x, k^) = glgg 9 {x) 6» (1 - x) ^ Q 

f 1 


2=0 


+ 


:(l-x) 


~ x{l — x) _|_ Mf — X (1 — x)]^ 


The integration over k^ of the TMD yields the pion PD 


(5) 


q{x) = J (Pkp {x, k^) 

= J ^ {p\ {C) 4 (oo, C) 7^ -Co (oo, 0) Ip (0) \p) , 

with A:+ = xP '^, = 0 and = 0. One gets explicitly 

q{x) = glqq 9 (x) 6» (1 - x) 

V.- I'ln I mlx{l-x) \ 

hiV [Mf-mlx{l-x)) [M2-TO2x(l-x)]j ■ 


( 6 ) 


(7) 


We stress that, since we are working in a field theoretical scheme, the right support of the 
distributions, 0 < x < 1, is not imposed and arises naturally. For the same reason, one can easily 
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proof that: i) the normalization is correct, i.e., f dx q{x) = 1 ; ii) J dx x q{x) = 0.5, i.e., the 
fraction of momentum carried by each quark is one half of the total momentum. Since at this level 
there are no sea quarks, this is the expected correct result. 

2.2 Boer-Mulders function 

The BM function is defined as 


i±u{d) 

'd,7r 


(x, — 


X 


^2k^J (2^)3 

(oo,^“) icr *+75 i(l±CT 3 ) Co (oo,0) V’(O) \p) 


( 8 ) 


with the same conventions used in Eq. O- To fix the ideas, as previously done for the unpolarized 
TMD, we will consider the BM function for a u quark in a 7 r+. 

At zero order in gs, vanishes, due to the T-odd character of the BM function. In order 
to have a non-zero value of the BM function, we expand the gauge link C^^ (oo,^“), Eq. @ in 
powers of gs, up to the first order, as it has been done for phenomenological model estimates of 
T-odd parton distributions, for the nucleon (see, e.g.. Refs. [IZl m [Ml [501E]) and, recently, for 
the pion m- For the gluonic field, we use its definition in terms of the source 

iv) = J d'^y D+’'(rj-y) gs^iy) > (9) 


where (x — y) is the gluon propagator. After some calculation, we arrive at 


loo, 


X 


-1 

-\- is 


(27r) 

J d'^y y 7"^ ^ (?/) 


t+ — ie 


( 10 ) 


At the first order in g1 one gets 




Ix/rp 

d'^t 


pitT^T p-it+C 


-1 


^ ( 27 r)"‘"" t+ - ie + ie 

J d^y (y) y 7 "^ '0 (y) * cr* + 75 
A“ r dH 1 


( 11 ) 


i + 

ter J5 


-1 


2 J ( 27 r)'‘ t+ - ie + ie 
A straightforward calculation leads to 


d'^y -tpiy) tlj{y)\ tlj{0)\p) . 


rJ- (t i, 2 \ _ 13 j j 2 2 

^T) — ^ 2 ysyiTqq 


dH dr 


(27r)^ A+ — ie -I- ie 


X {-Tr [Sf {r + ) 75 t- Sf [r + A+f) 

X S'p, (r 4 -f) 75 r+^F (r-f) ^7+] 

-I- Tr [5f (r-^) 75 T-5 'f (x-b-l) ^io-*+75 

X 5 'f (r-t-bf) 75 t+^f (r-t-f) ^7+]} , 


( 12 ) 


where r = ^(x — j . The traces in the equation above, in the order they appear, 

correspond to the diagrams in the left and right panels of Fig. [21 respectively. In principle, the BM 
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Figure 2: The diagrams describing the two traces in Eq. ([T2|). 


function could have contributions also from the sigma term (the one reported in Fig. [U right panel, 
in the unpolarized case). The direct calculation shows anyway that these contributions vanish. 
After a lengthy calculation, and including the Pauli Villars renormalization, we obtain 




X In 


fcy + Mf — x{l — x) 
Mf — toJ x{1 — x) 


where we have introduced the strong coupling constant = |^. 
The kx integral can be calculated, providing 


K{x) = J (fkTh^.^{x,kx) =--^^mm^asglgge{x) e{l-x) 


X E' 

i=0 

which, integrated over x, yields 


* Mf — x{\ — x) 


2 4 arccsc 


/ /" 1 n ai c 

dx J d?kT [x, kl) = y] c* 


(2^) 

\ rn„ J 


[Mf-mD 


In the ^ 0 limit, we get 

2 


(13) 


(14) 


(15) 


(16) 


At variance with the fi^nix, kf^) case, in which f dx dkf, /i_^(a;, k^) = 1 is a consequence of charge 
conservation, the quantity (1161) is not in general related to any physical observable. However we 
note that, in the present NJL framework and in the chiral limit, one has gf^^ X)i=o 
i.e., the right hand side of Eq. (ITB)) can be related to r^r, the charge radius of the pion. 


3 Discussion and comparison with lattice data 

3.1 Unpolarized TMD 

To have a pictorial representation of the global x and kx dependencies, 3D-plots are shown in Fig. 
[21 for niTr = 140 MeV (left panel) and m.^. = 518 MeV (right panel). In the left panel, it can 
be seen that the unpolarized TMD varies slowly with x. This is easily understood looking at Eq. 

where x dependent terms always appear multiplied by In the right panel of Fig. [3] it is 
clearly seen that, by taking a heavy pion with m,r = 518 MeV, a value which will be useful later 
for the comparison with lattice data, the x dependence becomes much more pronounced. In the 
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Figure 3: The unpolarized TMD /i^jr as a function of kx and x, with = 140 MeV (left 

panel) and ttItt = 518 MeV (right panel). 





Figure 4: The pion PD, Eq. ([7|), after LO (solid line) and NLO (dashed line) QCD evolution to 
Q = 2 GeV (picture taken from Ref. [23]). Data are from Ref. [M]. 


latter situation, our results agree qualitatively with the findings of Ref. [mile], where different 
constituent quark models have been used to evaluate the unpolarized TMDs. This fact can be 
understood thinking that, in the present NJL approach, the chiral limit is naturally included, 
at variance with a constituent quark scenario, where chiral symmetry is explicitly broken. As a 
consequence, the x dependence of our results with a pion mass of = 518 MeV is closer to that 
obtained within constituent quark models, with respect to what is obtained in our approach using 
the physical pion mass. 

The rather flat x dependence, obtained using = 140 MeV, is not a drawback of the model. 
Hadron models, like the NJL model, must be regarded as a realization of QCD at a very low . 
Evolution will change the x dependence in an important way. In fact, starting from Eq 0, in Refs. 
[32] and |23l |33] , a very good description of the data of the pion parton distribution at Q = 2 GeV 
[34] is obtained, as one can see in Fig. IH taken from [23]. For later convenience, it is useful to 
report that the LO parameters of the QCD evolution used in Refs. [32] and [23l|33] predict as(2 
GeV) = 0.32 and as(2 GeV) = 0.29, respectively. These values are in good agreement with Og 
measured in correspondence of the mass of the r lepton, as{mr = 1.777 GeV) = 0.327to;oi6 [55] . 

Concerning the relation between mass and x dependence, it is also interesting to observe that, 
in the chiral limit, the NJL model predicts an absolutely flat parton distribution, q^ix) = 1, and 
distribution amplitude, ipxi^) ~ 1- Nevertheless, the different regime of evolution (DGLAP for the 
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Figure 5: The unpolarized TMD (x, k^), Eq. ([5|), at x = 0.5, for the physical pion, as a function 
of kT (full). The dashed curve does not include the counter terms coming from the regularization 
procedure. 


first quantity and ERBL for the second one) produces very different x dependencies at higher 

m- 

In Fig. [5l the kx dependence is shown, having fixed x = 0.5. The result without the con¬ 
tributions of the counter terms originated by the regularization procedure is also reported. It is 
worth stressing that, in our approach, the kx dependence is automatically generated by the NJL 
dynamics. This is an important feature of our results, not found in other approaches. In facts, for 
example, the two different kx dependencies of the unpolarized TMD shown in Ref. m are dictated 
by two different forms adopted for a regulator function appearing in the pion Bethe-Salpeter am¬ 
plitude. In a similar fashion, in the Light-Front scenario of Ref. [16] , the obtained IcT-dependence 
is determined by the gaussian form assumed in the pion light-cone wave function, following Refs. 
IMlISZj- In our case, the kx dependence is not imposed using an educated guess. It is therefore rel¬ 
evant to report that our prediction has the following asymptotic kx behavior, as it can be obtained 
from Fq. 


/l, 7 r (x, krp^ 


fc|,^-oo 27r3fc6 


(17) 


We reiterate that this is just a consequence of the NJL model with Pauli-Villars regularization. To 
this respect. Fig. [S] points out the importance of the regularization procedure. In facts, without 
the counter terms, which suppress the high kx region, the TMD would not be integrable in the 
variable Actually, as it has been shown above, the integration over of the TMD yields the 
pion PD. 


3.2 Boer-Mulders function 

Numerical results of the evaluation of the BM function, Eq. (d, divided by the strong coupling 
constant, are reported in Fig. 16] in a 3D-plot, providing a pictorial representation of the global x 
and kx dependencies. 

As it happens for the unpolarized TMD, the BM TMD varies slowly with x when the physical 
pion mass, mjr = 140 MeV, is used in our calculation. This is again easily understood looking at 
Eq. (I13L where x dependent parts always appear multiplied by m^. In the right panel of Fig. [6l 
it is shown that, by taking tUtt = 518 MeV, the x dependence becomes much more relevant. 

As for the unpolarized TMD, the obtained kx behavior is a genuine result of the NJL dynamics 
with Pauli-Villars regularization. We report our prediction for the asymptotic kx behavior of the 
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Figure 6: The Boer-Mulders TMD, (x, , divided by the strong coupling constant a*, evaluated 

using = 140 MeV (left panel) and = 518 MeV (right panel). 


BM TMD, as it can be obtained from Eq. (1131) : 


^ ,2^ , n, ^ - mlx{l - x) 

hi,^[x,k-r) — -^ e{x)e{l-x) - > c,log— - — -- 

V 27 r 3 fc^ ^ m^-m^x{l-x) 


k^—¥oo 


+ O (^-^log{kT)^ 


(18) 


3.3 Comparison with lattice data 


In the following, we compare our results with lattice measurements. In facts, very recently, a lattice 
calculation has been performed [5] , focused on a TMD observable related to the Boer-Mulders effect 
in a pion. 

The quantity which has been addressed is a ratio, defined in an appropriate way to be safely 
evaluated on the lattice. It is the so called “generalized Boer-Mulders shift”, given by the following 
expression 


{ky)uT(bT) 




(19) 


where and are x-moments of generic Fourier-transformed TMDs: 


f[i](o) 


/H(-)(6^) = j^^dxx^-^ j f{x,kl) (20) 

= (55^ y^xx^-iydlfeTlIfeTl (^|^yj„(|6T||A:T|)/(x,fc^), 

with Jn denoting the Bessel functions of the first kind. 

One should notice that the &t —t 0 limit of the generalized Boer-Mulders shift, Eq. dUD, 
formally corresponds to fc|.-moments of TMDs, 

/MO)( 0 ) = J'^dxx^-^ J d^kr f{x,kl) ( 21 ) 

In the 6 t —>■ 0 limit, the generalized Boer-Mulders shift reduces therefore to the “Boer-Mulders 
shift”, 

(M^T(0)=m.-Tjj^, (22) 

which has the meaning of the average transverse momentum in y-direction of quarks polarized in 
the transverse (“T”) x-direction, in an unpolarized (“t/”) pion, normalized to the corresponding 
number of valence quarks. 
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Figure 7: The moment of the unpolarized pion TMD, calculated assuming the physical 

pion mass, m-,^ = 140 MeV (full line), or m,r = 518 MeV, the value used in the lattice calculation in 
Ref. [8] (dashed line). 



Figure 8: The same as in Fig. [Tj for m-,^ times the moment of the Boer-Mulders pion 

TMD, divided by a<j- 
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Figure 9: The generalized Boer-Mulders shift, Eq. (jl9p . as a function of The orange circles 
are obtained using the present NJL approach with m-,^ = 140 MeV; the red stars are obtained using 
mjT = 518 MeV, the value at which the lattice measurement has been performed. Two different sets of 
lattice data (black and blue boxes), with their RMS deviation, obtained in Ref. [8] using independent 
fits, are shown for comparison. 


It should be noted, however, that the fc^-moments of TMDs CD appearing in (1221) are in 
general divergent at large Ict and thus not well-defined without an additional regularization. In the 
generalized quantity, dUD, a finite bx effectively acts as a regulator through the associated Bessel 
weighting, cf. (l20l) . For these reasons, in Ref. [8], lattice QCD data have been obtained for the 
generalized Boer-Mulders shift (IT^ . at finite bx, using a pion mass of 518 MeV. 

In the following, we compare these lattice results with the outcome of our approach. The 
numerator and denominator of the generalized BM shift, Eq. (ITOl) . i.e., the moments 
and are shown in Figs. [7]and[8l respectively, as a function of bx, for = 140 

MeV (full line) and = 518 MeV (dashed line). In Fig. [Dwe observe that, by increasing the pion 
mass, f[^^^^\bx) is shifted towards higher values of bx- This is consistent with the fact that the 
e.m. radius of the pion is smaller for = 140 MeV than for = 518 MeV (see the Appendix 
for the actual values). In Fig. [8] we compensated an overall factor m~^ present in 
multiplying the latter quantity by the corresponding pion mass. We observe the same behavior, in 
relation with the variation of the mass, as in Fig. [T] Nevertheless, it is difficult to give any simple 
intuitive explanation because here we are dealing with a two body operator, as it can be seen from 
Fig. [2] or Eq. (HD- 

In Ref. [8], lattice data for the generalized Boer-Mulders shift have been presented for three 
different values of bx, at the momentum scale Q = 2 GeV, given in [38]. Our model results have 
therefore to be evolved to this scale, for a proper comparison with the lattice calculation. Unfor¬ 
tunately, the complete QCD evolution of the moments and , involving both 

the X and kx dependencies, is not under control. In particular, the evolution of the kx dependence 
is basically unknown. To estimate the evolution of the x dependence of the denominator, we can 
approximate it with the behavior of the corresponding kx integrated quantity, the first moment of 
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|6t| 

fm 

Lattice 1 

GeV 

Lattice 2 

GeV 

NJL = 140 MeV 

GeV 

NJL = 518 MeV 

GeV 

0.27 

-0.138(28) 

-0.133(19) 

-0.104 

-0.087 

0.34 

-0.128(29) 

-0.121(16) 

-0.093 

-0.076 

0.36 

-0.145(25) 

-0.148(15) 

-0.090 

-0.074 


Table 1: The generalized Boer-Mulders shift, Eq. (I19p . as a function of hx- The first column corre¬ 
sponds to the three values of hx analyzed in Ref. [8]. The second and third columns contain the two 
different sets of lattice data, with their RMS deviation, obtained in Ref. [5] using independent fits. 
In the fourth and fifth columns the results obtained using the present NJL approach are given, with 
= 140 MeV and = 518 MeV, the value at which the lattice measurement has been performed. 

the PD, which does not evolve in x. For the numerator, being proportional to as(Q'^) 

through the BM function, for a first estimate one can assume that its evolution is basically governed 
by that of as{Q‘^). It is therefore important to fix properly the value of as in evaluating the model 
prediction for the generalized Boer-Mulders shift at 2 GeV. Following the discussion on the fixing of 
the LO evolution parameters in NJL calculations of parton distributions, reported in the previous 
section, we use as = 0.31. 

In Fig. [9] and Tab. [T] our results are compared with the lattice data, evaluated according to 
two different fits providing consistent results [6]. We obtain a reasonably good agreement. It must 
be emphasized that our calculation has been performed in the NJL model without introducing any 
new parameter. We observe that the generalized Boer-Mulders shift varies slowly when we go from 
ttItt = 140 MeV to = 518 MeV. 

The main uncertainty in our calculation comes from the poorly known QCD evolution of the 
moments of the TMDs, entering the definition of the generalized Boer-Mulders shift. Summarizing 
our approximated evolution scheme, the x evolution of the denominator has been neglected thinking 
to the behavior of the corresponding PD, the one of the numerator has been assumed to be governed 
by that of as only, and the kx evolution has been neglected overall. 

4 Conclusions 

We have considered the well-established NJL model, without any additional parameter, for the 
study of the two leading twist pion TMDs, the unpolarized, and the Boer-Mulders 

one, hj-.^{x, k^). We were motivated by the success of this model in reproducing pion observables, 
such as the parton distribution and the pion gamma transition form factor, and by the aim of 
reproducing recent lattice results [8]. Since in the latter calculation a value of 518 MeV has been 
used for the pion mass, we present our results for = 140 MeV and m,r = 518 MeV. 

We have studied the kx dependence of and hi^^{x,k'^). In both cases, this depen¬ 

dence is automatically generated by the NJL dynamics. The obtained kx asymptotic behavior 
of these two quantities, at the momentum scale of the model, Qoi is found to be and 
respectively. Nevertheless, QCD evolution to higher scales could modify these trends. 

We observe a soft dependence on x in both TMDs at Qq. This can be easily understood 
observing that, in the final expressions of the TMDs, the x-dependent part is always multiplied by 
ttitt- Our experience with the parton distribution and the distribution amplitude of the pion is that 
this X dependence provides remarkably good results after QCD evolution. When the = 518 
MeV case is considered, we get a stronger x dependence, approaching results obtained in models 
built with constituent quarks. 

Finally, we have studied the generalized Boer-Mulders shift, which has been recently calculated. 

The agreement we obtain with these lattice data is rather good, qualitatively and quantitatively. 

Our results show a weak dependence on the mass of the pion. 

The main theoretical uncertainty in our calculation comes from the approximated QCD evolu¬ 
tion we have performed. A more conclusive comparison would require therefore further lattice data 
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and the implementation of the correct evolution of the TMDs moments appearing in the calculation. 
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A The NJL model and regularization scheme 


The Lagrangian density in the two-flavor version of the NJL model with electromagnetic (e.m.) 
coupling is 

C = '4){i p - mo) tp + g 

with D^ = d^+ie A^. 

The NJL model is a non-renormalizable field theory and a cut-off procedure must be defined. 
We have used the Pauli-Villars regularization in order to render the occurring integrals finite. This 
means that, for any integral, we make the replacement 


+ ('0Tf75 i’Y 


Pq 


f{g\m) 


Pq 


J2c,f{q;M,) 

i=o 


with Mj = rp -I-j A^, Co = C2 = 1 and ci = —2. Following ref. m we determine the regularization 
parameter A and m by calculating the pion decay constant and the quark condensate in the chiral 
limit, via 


2 3 m \ 

j=o 


Mf 

log^ , 


(uu) 


2 


-T-7 > . log 


4jj-2 

j=o 


2 


with mo fixing the pion mass. 

With the conventional values {uu) = —(0.250GeV)^, /tt = 0.0924GeV and m^ = 0.140GeV, 
we get m = 0.238 GeV, A = 0.860 GeV and mo = 5.4 MeV. For the pion-quarks coupling constant 
we get = 6.279. The electromagnetic pion radius turns out to be rj = 0.31 fm^ (experimental 
value rj = 0.44 fm^). 

For a proper comparison with lattice data, we have applied the same model to a massive pion, 
with mTT = 0.518 GeV. We have not changed the value of A; for m we have taken m = 0.300 GeV. 
In this way we have {uu) = —(0.263 GeV)^, /tt = 0.0997 GeV, mo = 83 MeV and = 3.667. The 
e.m. pion radius is = 0.38 fm^ . 
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